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PERFORMANCE  ANALYSIS  OF  A SPLIT-BEAM 
AUTOMATIC  TARGET  FOLLOWER  IN 
ADVERSE  ACOUSTIC  ENVIRONMENTS 


INTRODUCTION 

Present  performance  analyses  of  broadband  split-beam  bearing  track- 
ers can  be  grouped  into  two  categories: 

1.  the  open-loop  analysis  of  trackers  searching  for  a stationary 
target  in  a relatively  general  noise  field1’4  and 

2.  the  closed-loop  analysis  of  trackers  searching  for  a dynamic 
target  in  a simple,  specific  noise  field.5"7 

The  results  of  reports  referenced  under  category  1 generally  yield 
performance  bounds,  that  is,  values  of  bias  and  random  errors  under  a 
specific  set  of  acoustic  conditions  that  can  only  be  approached  by  any 
practical  automatic  target  follower  (ATF). 

On  the  other  hand,  those  referenced  under  category  2 yield  design 
criteria  on  the  feedback  loop.  By  linearizing  the  bearing  deviation 
indicator  (BDI)  about  the  expected  operating  point  and  using  conven- 
tional control  theory,  a conceptually  simple  model  is  obtained  for  a 
very  complex  control  problem.  Unfortunately,  the  model  is  strictly 
applicable  only  to  a single  moving  target  in  an  uncorrelated  ambient 
noise  background. 

In  contrast,  this  report  presents  a closed-loop  analysis  of  an  ATF 
that  is  searching  for  one  or  more  stationary  targets  in  a relatively 
general  noise  field.  The  analytical  technique  differs  from  that  of 
category  2 in  one  important  respect.  Whereas  the  controlled  delay  (e) 
in  the  BDI  shown  in  figure  1 is  initially  (before  linearization)  as- 
sumed to  be  deterministic,  e is  considered  to  be  a random  variable 
in  this  report. 

As  a result  of  this  more  general  approach,  two  coupled,  transcen- 
dental equations  involving  the  mean  and  variance  of  e are  obtained. 

A simultaneous  graphic  solution  of  these  equations  will  then  yield 
specific  values  of  the  bias  and  random  errors  for  any  given  adverse 
acoustic  environment. 
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In  order  to  focus  on  the  ATF  analysis,  only  a two-element  array 
is  postulated.  However,  the  derivation  described  in  appendix  A may 
readily  be  modified  to  include  any  array  geometry. 


DISCUSSION— BDI  WITH  RANDOM  TIME  DELAY 

The  BDI  considered  in  this  report  is  the  analoq  anticorrelator 
shown  in  figure  1,  where  the  inputs  Wi(t)  and  w2(t)  are  assumed  to 
be  sample  functions  of  a stationary  zero  mean  Gaussian  random  process 
Filters  h(t)  and  y(t)  are  real;  additionally,  y(t)  must  be  odd*  so 
that  the  device  is,  in  fact,  a deviation  indicator.  As  indicated  in 
figure  1 , a 90°  phase  shifter  or  a differentiator  satisfies  this  lat- 
ter requirement.  Finally,  the  random  time  delay  (e)  is  assumed  to 
be  uncorrelated  from  wx(t)  and  w2(t). 

In  appendix  A,  the  mean  (m  ) and  the  variance  (o^)  of  the  BDI 
output  u are  derived  in  terms  of  the  characteristic  function  of  an 
arbitrary  random  time  delay  (e).  It  is  further  shown  that  if  e is 
Gaussian  with  mean  m and  variance  a2,  then 


m, 


iuK<,2)  . /df  e-V2(2.f») 

2(m.o2)  = E -2n 


n=l,2,...  nT 


Auf  + f GNN(f,m,a2)df, 
8mn  / 


(1) 

(2) 


where 


G^(f ) = 


♦,(f)  © *2(f)  +/dx  4>1 2(x)  4>-j g (x  - f)ei27rm(2x  ' f) 


• e~  1 / 2[ 2tts ( 2x  - f)]‘ 


(3) 


*A1 ternatively,  the  requirement  on  y(t)  is  that  its  transform  (Y(f)) 
must  be  purely  imaginary  and  odd.  Y(f)  = Y*(-f)  = -Y(-f). 
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♦i(0  = |H(f)r  G (f)  m = T 

i w1w] 

*2(f)  « |H(f)|2  | Y(f ) | 2 (f)  o2  = e2  - e 

*12(f)  » |H(f)|2  Y*(f)  (f), 

where 

® = denotes  convolution, 

G (f)  = auto-spectral  density  of  input  w, (t), 
wlwl  1 

G..  w (^)  = auto-spectral  density  of  input  w~(t), 
w 2*'  ^ ^ 

Gw  w (f)  = cross-spectral  density  of  inputs  w,(t)  and  w,(t),  and 
w ^ w 2 i c 

H ( f ) and  Y(f)  = Fourier  transforms  of  filter  h(t)  and  y(t),  respectively. 


As  noted  in  appendix  A,  GNN(f)  represents  the  average  (averaged 
over  e)  auto-spectral  densityMof  a random  variable  (N),  which,  at  a 
given  e,  is  the  deviation  of  u from  its  mean  (u) . It  is  a generalized 
form  of  the  equivalent  noise  source  or  disturbing  perturbation  used 
by  Lindgren6  and  others  to  account  for  the  finite  averaging  time  of 
the  closed-loop  filter. 

In  order  to  appreciate  the  effect  of  the  randomness  of  e on  the 
output  (u),  consider  the  mean  and  variance  of  u when  the  "jitter"  is 
zero,  that  is,  when  e is  deterministic  with  value  m and  a = 0.  Under 
this  condition. 


mu  =/df  *12(f)  e12rfm  (4) 
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Notice  that  the  infinite  series  in  equation  (2)  vanishes  when  o = 0. 
The  reader  may  confirm  that  these  are  the  usual  expressions  found  in 
the  literature  for  the  analog  BDI  shown  in  figure  1. 

As  the  jitter  increases v(o  f 0),  mu  will  decrease  from  its  maximum 
value  given  in  equation  (4),  whereas  ay2  will  generally  increase 
because  the  contribution  from  the  infinite  series  in  equation  (2)  is 
greater  than  the  decrease  in  G 

the  jitter  becomes  very  large,  . . , 

ever,  o2  approaches  some  limiting  value.  The  interpretation  of  this 
behavior  is  that  the  output  signal-to-noise  ratio  (SNR)  approaches 
zero. 


ATF  ANALYSIS 

The  complete  ATF  is  formed  by  "feeding"  the  BDI  output  (u( t) ) into 
an  appropriate  feedback  filter  (z(t) ) and  using  its  output  (v(t))  to 
control  the  delay  ( e ( t ) ) . Figure  2 shows  the  resultant  circuit  to 
be  analyzed.  The  dashed  line  represents  that  portion  of  the  feedback 
loop  that  will  be  closed  in  the  following  discussion. 

It  is  assumed  in  this  analysis  that  the  effective  integration 
time  of  the  closed  loop  is  much  longer  than  the  correlation  time  of 
the  inputs.  Thereby,  when  the  loop  is  closed,  v(t)  approaches  a sim- 
ilar (but  not  the  same)  Gaussian  distribution  as  that  assumed  for  e 
in  the  derivation  of  equations  (1)  and  (2). 

STEADY-STATE  MEAN  CALCULATION 

If  we  assume  that  the  loop  is  open,  the  expectation  (m  ) of  the 
random  variable  (v)  is  given  by 


NN(f)  (see  equation  3).  Finally,  as 
m aDDroaches  zero  for  all  inouts:  how- 


my  = Z(0)mu,  (6) 

where  Z(0)  is  the  dc  gain  of  the  feedback  filter  whose  units  are  such 
that  v has  units  of  time  (sec). 

The  loop  is  now  closed  in  the  usual  manner  by  stipulating  that, 
in  the  steady  state,  the  mean  of  e be  equal  to  the  mean  of  v.  Using 
equation  (1)  and  denoting  the  steady-state  values  of  m and  a as  mQ 
and  o , respectively,  yields 
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I 


m 

o 


m = Z(0)m 
v ' ' u 


Z(0)  J df 


$12(f)e 


i2iTfm 


■0e-l/2(2nfo0)< 


(7) 


which  is  a transcendental  equation  in  m whose  solutions  represent 
possible  steady-state  values  of  m for  any  given  value  of  cqz. 

Typically,  the  feedback  filter  ( Z ( f ) ) is  at  least  of  the  first 
order,  so  that  Z (0 ) tends  to  infinity.  Equation  (7)  then  reduces  to 

"u-/df*  12(f)e,2*f"#  .-1/2(2,f'.)2  (8) 

Thus,  for  a given  oQ,  the  possible  steady-state  values  of  m are 
given  by  those  values  for  which  the  mean  BDI  output  (error  function) 
is  zero,  e.g.,  by  the  zero  crossing  of  the  error  function. 

Notice  that  the  derivation  of  equation  (7)  or  (8)  has  not  required 
any  linearization.  Furthermore,  these  expressions  are  valid  for  any 
stationary  inputs  (wj(t)  and  w2(t))  provided  the  effective  integration 
time  of  the  loop  filter  is  large  with  respect  to  the  correlation  time 
of  the  inputs. 


STEADY-STATE  RANDOM  ERROR  CALCULATION 

The  small  steady-state  random  deviations  of  e may  be  obtained  by 
the  usual  practice  of  linearizing  the  BDI  and  using  control  theory  to 
derive  the  closed-loop  transfer  function.  For  this  purpose,  assume 
that  the  ATF  has  achieved  a steady-state  condition,  so  that  the  fluc- 
tuations in  the  BDI  output  are  only  due  to  the  jitter  of  e.  Denoting 
the  steady-state  values  of  m and  o by  m and  o , respectively,  we 
see  in  appendix  B that  the  variance  of  the  BDI°output  about  this 
operating  point  is  given  approximately  by 


°u(Va2)  ^ 0 KBDI^mo’ao^  +/WfVo>-. 


(9) 


where 


(m  a2) 

' o o' 


(10) 
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Now,  in  the  steady  state,  the  mean  of  the  BD I output  is  zero  and  its 
variance  is  given  by  equation  (9).  A linear  system  that  matches  these 
requirements  is  shown  in  figure  3 if  the  mean  of  the  equivalent  noise 
source  (N1)  is  zero  and  if  its  spectral  density  is  (f,  m,  a2). 

By  convention,  KBQj(mQ,ap)  is  denoted  as  the  "BDI  gain."  With 
reference  to  eauation  (lOj,  the  gain  is  equal  to  the  total  derivative 
(variation)  of  the  error  function  and,  therefore,  represents  a more 
general  form  of  the  conventional  BDI  gain  expression.  Further  compari- 
son shows  that  when  c is  small,  the  first  term  in  (10)  dominates, 
and  Kbq i reduces  (forzero  bias  error)  to  the  expression  obtained  by 
Pinkos  and  Develet.0 

It  should  be  noted  that  the  linearized  model  of  the  BDI,  as  shown 
in  figure  3,  is  applicable  only  at  the  steady-state  operating  point, 
and  its  usefulness  is  limited  to  determining  the  steady-state  jitter 

(«.>• 

Accordingly,  after  passing  the  BDI  output  (u)  through  the  feedback 
filter  (z),  the  ATF  loop  may  be  closed,  as  shown  in  figure  4(a).  Via 
conventional  control  theory,  the  equivalent  closed-loop  filter  (A(f)), 
which  relates  the  input  (N'(t))  to  the  output  (e(t)),  is  given  by 

ft(-C)  . z(f)  . (11) 

‘ ’ ' ’ * Sji  Z(f) 

Finally,  with  reference  to  figure  4(b),  and  remembering  that  N‘  = 0, 
the  variance  of  z becomes 

2 

°o2  ■ /df  S'N-Cf)  */df  <Wf>  1 ■*  021 

Under  the  initial  assumption  of  long  integration  times,  the  band- 
width of  ) A ) 2 will  be  much  smaller  than  that  of  GN • N « » so  that,  to 
a good  approximation,  the  variance  is  given  by 


2 _ GNN^0,mo,(V 

* I’Wv’M2 


wf6 


' + kbdi(V°o>  z(f> 


r 
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where,  from  equation  (3), 


GNN(0,mo’°o)  = / df  *2(f)  + *?2(f)  el47Tfm°  e'1/2(47Tfoo)2]- 

(14) 

Substituting  for  fi, , and  making  use  of  the  restrictions 

imposed  on  Y(f)  yields^  'c 

GNM<°-,V’o>  = /df  IW)!2  [G„iW](f)  G^f) 

-Gw,w2(f)e(4’f,no  e-l/2(4,fo0)2J.  (15) 

Typically,  one  defines  the  integral  in  equation  (13)  as  the  equiv- 
alent closed-loop  bandwidth  (Afg);  that  is, 

, ? f KBDI(mo’°o^  Z(f)  2 

Af  (m  ,aZ)  J J-D-1  .0  o df.  (16) 

6 ° ° J 1 + KBDI(Vao>  Z<f> 

Thereby,  equation  (13)  becomes 


2 ~ Sn^’V'V 


°o  ” 


KBDI^mo,0oM 


AfJrn  ,0 
e o o 


Notice  that  equation  (17)  is  also  a transcendental  equation  in 
a whose  solutions,  for  a given  value  of  m0,  are  possible  steady-state 
values  of  the  random  error.  Finally,  the  simul taneous  solutions  to 
equations  (8)  and  (17)  represent  possible  steady-state  values  for  the 
mean  and  variance  of  the  random  time  delay  (e).  Such  an  analysis 
for  a typical  ATF  and  for  four  adverse  acoustic  environments  is  pre- 
sented in  the  next  section. 


iOnai 
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TYPICAL  RESULTS 

In  this  section,  the  steady-state  bias  and  random  errors  are  cal- 
culated for  a typical  ATF  that  is  operating  in  four  adverse  acoustic 
environments: 

Case  I — single  target  in  uncorrelated  noise 

Case  II  — single  target  in  partially  correlated  noise 

Case  III  --  single  target  in  uncorrelated  noise  plus  interference 

Case  IV  — single  target  with  multipath  in  uncorrelated  noise 

In  all  these  acoustic  environments,  the  signal,  noise,  and  interfer- 
ence spectral  densities  are  assumed  to  be  uniform;  the  noise  power 
is  assumed  to  be  equal  at  each  input;  and  the  direct-path  signal 
spectral  density  is  assumed  to  be  unity. 

With  reference  to  figure  5,  the  general  ATF  inputs  Wj ( t)  and  w2(t) 
for  the  four  environments  are 

Wj(t)  = S(t  + t0$1)  + 6S(t  + tQS2)  + I(t  + tQI)  + n^t) 

w2(t)  = s(t  + tQS1  + t$1)  + es(t  + tQS2  + t$2) 

+ I(t  + tQI  + tj)  + n2(t),  (18) 

where 

S 

65 

I 

n-j  i n2 

^Sl 
l0S2 


= desired  signal  (direct  path) 

= multipath  signal 
= interfering  signal 
= noise 

= bulk  time  delay  of  the  direct  signal 
= bulk  time  delay  of  the  multipath  signal 
= bulk  time  delay  of  the  interfering  signal 
= jr  sines1  - interelement  delay  of  the  direct  signal 
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e 

e 

e 


51 

52 

I 


6 


jr  sines2  - interelement  delay  of  the  multipath  signal 
£ sinGj  - interelement  delay  of  the  Interfering  signal 
sensor  spacing 
speed  of  sound  in  water 

arrival  angle  of  the  direct  signal  with  respect  to  broadside 

arrival  angle  of  the  multipath  signal  with  respect  to  broad- 
side 

arrival  angle  of  the  interfering  signal  with  respect  to 
broadside 

amplitude  ratio  of  the  multipath  signal  with  respect  to 
the  direct  signal 


Accordingly,  the  auto-  and  cross-spectral  densities  of  wj(t)  and 
w2(t)  are 


Gw,w2<f>  ' 1 + **  + 26  2"  fto  + GII  * Gn„  <’9> 

Gw2w2(f)  = 1 + e2  + 23  cos  2fff (tQ  + tsl  - ts2)  + Gn  + Gnn 

(20) 

G (f)  = e"l2rftSl  + Be"l27lf^tS2  ' V + Be‘l2*f(tSl  + V 
wlw2 

+ B2e"l2,lftS2  + Gjj  e'l27rftI  + y(f)Gnn  (21) 


where  the  previous  assumptions  regarding  the  acoustic  fields  have 
been  incorporated  and  where 


' t0S2 
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G n (f) 
n-i  n„ 

y(0  = Yn  n (f)  = -M 

"ln2  Gnn 


noise  coherence. 


Notice  that  for  the  above  assumptions,  G j j and  G_n  may  be  inter- 
preted as  the  interference- to-signal  and  noise- to-signal  ratio, 
respectively. 

The  particular  ATF  configuration  considered  for  these  four  acoustic 
environment  cases  is  shown  in  figure  6.  Additional  particulars  are 
as  follows: 


I H ( f ) | 1 


2 = 1 402.65  < |f|  < 805.3 

jO  otherwise 


Y(f)  = i 2-rrf  (differentiator) 


ot  (1  + i27TfT/3)  (1  + i2nfT/18) 

! i2rrf ) (1  + i27rfT)  (1  + i2^fT/9)  (1  + i2^fT/90) 


where  T = 50  sec  and  a is  a proportionality  constant  having  units  of 
sec/volt2  whose  value  must  be  chosen  in  accordance  with  dynamic  follow- 
ing considerations . 

Substituting  Z(f)  into  equation  (16)  and  performing  the  integra- 
tion yields 


Afe  = 


(«KnnTT)  Tl  + 0.11(aKDnTT)  + 0.001 (aKDnTT)2 


1 + 0.30(aKBDIT)  + 0.0067(aKBDIT) 


BDI_  I,  (22) 


which,  to  a good  approximation,  is 
(aKBDIT)0"7 

Afe  « — ^ for  1 < UKbdiT)  < 1000. 


(aKBDIT)( 
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The  restriction  on  a and,  therefore,  on  the  closed-loop  bandwidth 
is  determined  by  the  maximum  tolerable  "lag  distance."  With  reference 
to  equation  (16),  the  velocity  error  constant  (K  ) is  defined  as 


kbdi“ 


(24) 


As  shown  in  a Raytheon  study,5  a maximum  allowable  lag  distance 
of  100  yd  (91.44m)  for  a target  that  is  moving  at  a constant  speed  of 
40  knots  requires  that  K > 0.225.  Therefore,  the  value  of  a must  be 
chosen  so  that,  within  the  operating  range  of  the  ATF, 


a > 


0.225 

^di 


or 


(aKBDIT) 


> 11.25. 


(25) 


In  order  to  demonstrate  the  procedure  required  to  calibrate  the 
ATF  (chose  a)  and  to  solve  for  the  steady-state  values  of  m and  a, 
case  I will  be  worked  in  detail.  For  the  other  examples,  only  the 
final  results  will  be  presented. 


CASE  I— SINGLE  TARGET  IN  UNCORRELATED  NOISE 
The  conditions  for  case  I are 
e = 0,  I = 0,  y = 0, 
so  that  equations  (19)- (21)  reduce  to 

<L  „ (0  ■ G...  „ <f)  • 1 + G. 


"l"l 


W2W2 


nn 


(26) 


Gw  w (f)  = e"l2irfts1 
wlw2 


Substituting  equation  (27)  and  the  expressions  for  Y(f)  and  |H(f)| 
into  equation  (8)  yields 


(27) 

2 


2 

nu  = * / “ sin  [u(m0  - ts1)]  e 


-1/2(<O0  r da)  = 0,  (28) 


“1 
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where 


u = 2nf , uj i = 2n(402 . 65)  rad/sec,  wg  = 2tt(805.3)  rad/sec. 


The  error  function  (mu)  is  plotted  in  figure  7 as  a function  of 
the  bias  error  (b  = m - t,., ) for  zero  jitter  and  for  a0  = 0.25 
msec.  In  accordance  with  equation  (28),  possible  solutions  for  b are 
given  by  the  zero-crossinqs  of  this  function.  The  locus  of  these 
solutions  (zero-crossings)  are  shown  as  dotted  lines  in  figure  8.  It 
will  be  demonstrated  shortly  that  some  of  these  solutions  are  unstable 
while  others  cannot  exist  in  the  physical  ATF  due  to  the  incorrect 
sign  of  the  error  function. 

The  gain  of  the  BDI  is  obtained  from  equation  (10)  and  the  left 
side  of  equation  (28) : 


(29) 

which  is  plotted  in  figures  9 and  10  as  a function  of  bias  error  and 
random  error,  respectively. 


Notice  in  figure  9 that,  for  fixed  jitter,  the  gain  of  the  BDI 
cycles  with  the  bias  error,  with  maxima  occurring  in  the  vicinity  of 
the  zero-crossings  of  the  error  function.  Also,  notice  that  the  high- 
est gain  is  always  obtained  at  zero  bias  error. 


The  important  feature  to  be  observed  in  figure  10  is  that,  for 
small  fixed  bias  errors,  the  BDI  gain  curves  seem  to  "invert"  about 
the  point  (T).  That  is',  for  a < 0.25  msec,  the  gain  decreases  as  b 
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increases;  whereas,  for  a > 0.25  msec,  the  gain  increases  as  b 
increases.  The  significance  of  this  phenomenon  will  become  apparent 
shortly. 

For  the  purpose  of  discussing  system  stability,  it  becomes  expedi- 
ent to  use  the  approximate  expression  (23)  for  Affi.  Substituting 
this  into  equation  (17)  and  rearranging  terms  yields 


1.3  2 

KBDI^mo’°o^  = 


UT) 

2T 


0.7 


(30) 


Also,  substituting  equation  (26)  and  the  expression  for  Y ( f ) and 
I H(f ) ( 2 into  equation  (15)  yields 


Sn^’V0^ 


= l/2dw  J[(1 


+ G ) 
nn' 


_e_1/2(2u)co)  cos  2w(mo  - tsl)J  • (31) 

For  low  SNR's,  the  auto-spectral  density  will  become  independent  of  m 

and  a , leaving  0 

o 3 


2 7u)2  Gnn 

G-  (0*V°o>  58  -kr11 


NN 


for  G >>  1 . 
nn 


(32) 


Hence,  for  low  SNR,  the  right  side  of  equation  (30)  will  be  a constant 
(horizontal  line)  proportional  to  G2n< 

The  graphic  solution  to  equation  (30)  is  shown  in  figure  11,  where 
the  left  side  of  equation  (30)  is  represented  by  solid  curves  and  the 
right  side  of  (30)  is  indicated,  for  some  arbitrary  G and  a,  by  the 
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dashed  horizontal  line.  In  accordance  with  (30),  the  intersections 
of  the  horizontal  line  with  the  solid  curves  represent  possible  steady- 
state  solutions  for  a.  As  the  SNR  decreases,  the  horizontal  line 
will  shift  upward  and  the  new  intersections  represent  possible  solu- 
tions at  the  higher  noise  level. 

The  question  of  system  instability  may  be  answered  by  proper  inter- 
pretation of  figure  11.  Consider  the  operating  point  at  Qj(o  =0.2 
msec,  b = 0.15  msec).  Suppose  that  the  SNR  is  fixed  and  that°there 
occurs  an  incremental  increase  in  b.  The  new  operating  point  will  be 
located  on  the  horizontal  line  and  to  the  right  of  Qx;  that  is,  an 
incremental  increase  in  the  bias  results  in  a operating  point  Qx' 
whose  jitter  is  higher  than  that  of  Qx.  Similar  considerations  at 
operating  point  Q2  cause  a shift  to  the  left;  that  is,  at  Q2  an 
increase  in  bias  yields  a decrease  in  jitter.  Clearly,  Q2  is  an 
unstable  operating  point.  One  may,  therefore,  establish  the  following 
stability  criterion: 

1.  An  operating  point  is  considered  to  be  stable  if,  at  a fixed 
SNR,  an  increase  in  b yields  an  increase  in  o. 

2.  An  operating  point  is  considered  to  be  unstable  if,  at  a fixed 
SNR,  an  increase  in  b yields  a decrease  in  a. 

An  operating  point  where  an  incremental  increase  in  the  bias  error 
results  in  no  change  in  o represents  a point  on  the  threshold  of 
stability/instability.  Point  T,  discussed  in  connection  with  figure 
10  and  also  indicated  in  figure  11,  is  one  such  operating  point.  In 
general,  the  coordinates  of  the  threshold  operating  points  are  given 
by  the  location  of  the  extremes  (maxima  and  minima)  of  figure  9. 

Figure  8 shows  the  locus  of  these  threshold  operating  points 
(solid  curves),  as  well  as  the  locus  of  the  zero-crossings  of  the  error 
function  (dashed  curves).  The  shaded  regions  of  this  figure  represent 
unstable  operating  regions,  in  accordance  with  the  previously  defined 
stability  criterion.  Also,  figure  8 shows  four  possible  stable  oper- 
ating modes  of  the  ATF  for  a single  target  in  uncorrelated  noise  at 
low  SNR's.  From  figure  7,  two  of  these  modes  (at  b»  0.8  msec  and 
b «2.25  msec)  may  be  discounted  since  the  sign  of  the  error  function 
is  such  as  to  drive  the  operating  point  away  from  these  zero-crossings. 
Alternatively,  one  may  state  that  the  actual  ATF  will  seek  those  zero- 
crossings  where  the  slope  of  the  error  function  is  positive.  As  a 
result,  the  first  two  operating  modes  of  the  actual  ATF  (for  case  I) 
are 
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b = 0;  0 < |o  | < 0.250  msec  (lowest  mode) 

and 

1.57  msec  < |b|  < 1.65  msec  0 < |oQ|  < 0.257  msec. 


Notice  that  the  smallest  total  rms  error  is  obtained  when  the  ATF  is 
operating  in  its  lowest  mode. 

As  to  which  mode  the  ATF  will  actually  assume  depends  on  the  SNR 
and  the  initial  bias  conditions  when  the  ATF  is  activated.  Thus,  it 
is  very  possible  that  at  a high  SNR  and  a large  initial  bias  error  the 
ATF  will  slip  into  the  higher  operating  mode.  However,  as  the  SNR 
decreases,  only  the  lowest  operating  mode  is  possible.  This  is  evi- 
dent from  figure  11  since,  as  the  horizontal  line  shifts  upward,  there 
will  be  no  intersection  with  the  curve  representing  b = 1 .£  msec. 
Alternatively,  since  the  BDI  gain  decreases  as  the  bias  increases  (see 
figure  9),  a point  will  be  reached  where  equation  (30)  can  no  longer 
be  satisfied  unless  the  ATF  slips  into  the  lowest  operating  mode. 
Continuing  with  this  reasoning,  as  the  SNR  decreases  further,  we  note 
that  a point  will  be  reached  where  this  lowest  mode  will  also  become 
unstable.  This  threshold  point  is  indicated  as  point  T in  figures  8, 
10,  and  11.  The  various  parameter  values  associated  with  this  signifi- 
cant condition  are 

b = bT  = 0 

aQ  = o-j.  = 0.25  msec 

kbdi  ■ kbdi|t  * °-717  * ,0'° 

G = G \ = <SNR>T  = 0.0143  ’ or  ("18.46  dB). 
nn  nn  |y 

Thus,  based  on  this  graphic  analysis  and  assuming  that  the  ATF  is 
operating  in  its  lowest  mode,  the  ATF  will  become  unstable  when  the 
SNR  approaches  -18.46  dB.  At  the  verge  of  threshold,  the  bias  error 
will  be  zero,  and  the  random  error  will  be  0.25  msec. 


Having  determined  the  final  threshold  operating  point,  we  can  now 
choose  a.  Using  equation  (25)  and  the  BDI  gain  at  the  threshold  yields 


a = 3.137  x 10-11  sec/volt^  . 


Theoretically,  this  value  of  a ensures  that  under  uncorrelated  noise 
conditions  the  ATF  will  be  able  to  follow  a 40  knot  moving  target 
over  its  entire  operating  range  with  a lag  error  not  exceeding 
100  yd  (91.44  m). 

Once  the  system  is  calibrated  (chosen  a),  it  is  possible  to  solve 
equation  (30)  or  (17)  for  any  other  SNR.  Indeed,  once  a has  been 
determined,  one  may  program  equations  (8)  and  (17)  on  a computer  and 
initiate  an  automatic  two-dimensional  search  for  the  steady-state 
solutions  of  the  bias  and  random  errors  for  any  arbitrary  acoustic 
conditions.  The  results  of  such  an  automated  tracking  program  for 
uncorrelated  noise  are  shown  in  figure  12  as  a function  of  the  SNR. 

Only  the  random  error  for  the  two  lowest  operating  modes  are  shown. 

The  bias  error  is,  of  course,  zero  for  the  lowest  mode  and  varies 
between  1.57  and  1.63  msec  for  the  2nd  order  mode  (see  figure  8). 
Inspection  of  figure  12  shows  that  the  threshold  of  the  lowest  order 
operating  mode  (thresholds)  is  at  SNR  = -18.6  dB  and  that  the  random 
error  at  this  point  is  0.28  msec.  The  reason  that  both  of  these  values 
differ  from  the  expected  -18.46  dB  and  0.25  msec,  respectively,  is  that 
the  automated  program  uses  the  exact  equivalent  bandwidth,  Afe  (equa- 
tion 22),  rather  than  approximation  (23),  which  was  used  in  the  previous 
graphic  analysis.  As  expected,  the  2nd  order  operating  mode  thresholds 
at  a higher  SNR  (-16.5  dB)  and  exhibits  a larger  total  rms  error. 

For  comparison  purposes,  figure  12  also  shows  the  Cramer-Rao  lower 
bound  (CRLB)  for  a two-element  analog  BDI  followed  by  a 20  sec  averager 
(open-loop  tracker).  The  absolute  values  of  the  random  errors  cannot 
be  compared  directly  because  one  curve  pertains  to  a closed-loop  ATF 
and  the  CRLB  pertains  to  an  open-loop  ATF.  What  is  significant,  how- 
ever, is  the  shape  of  the  two  curves.  The  figure  shows  that  o0  is 
proportional  to  the  theoretical  variation  over  a large  range  of  SNR: 


1 + 2 (SNR) 


L (SNR) 


T 


1/2 


At  a low  SNR,  the  ATF  exhibits  the  expected  threshold;  however,  at  a 
high  SNR,  the  model  gradually  becomes  invalid  due  to  the  assumption 
that  Gnn  (0,  m,a2)  is  independent  of  m and  o. 
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CASE  II— SINGLE  TARGET  IN  PARTIALLY  CORRELATED  NOISE 
The  conditions  for  case  II  are 
B = 0,  1=0, 

so  that  equations  (19)-(21)  reduce  to 


wlwl 


(f)  = G ,,  (f)  = 1 + G 

w^w^  nn 


G w (f)  = e'l27lftSl  + y(f)  Gnn" 

wlw2 


(33) 


These  spectral  densities  may  now  be  substituted  into  the  appro- 
priate expressions  in  the  manner  outlined  for  case  I.  When  the  pre- 
viously determined  value  for  a and  a direct  signal  time  delay  of 
t$i  = -0.2  msec  are  employed,  the  bias  and  random  errors  are  shown 
in  figure  13  as  a function  of  SNR,  with  the  noise  coherence  (y)  as 
a parameter.  For  this  case,  the  noise  coherence  is  assumed  to  be 
real  and  independent  of  frequency.*  Furthermore,  only  the  results  for 
the  lowest  operating  mode  are  shown. 

Notice,  when  comparing  figures  13  and  12,  that  as  the  noise 
becomes  partially  correlated,  the  bias  error  increases  substantially. 
As  the  SNR  decreases,  the  ATF  will  tend  to  track  the  noise,  which  is 
effectively  located  at  broadside.  Eventually,  the  bias  error  ap- 
proaches ( - tsi ) =0.2  msec.  Note  that  for  all  SNR's,  the  random  error 
is  only  a small  fraction  of  the  bias  error.  Although  not  shown  in 
figure  13,  the  ATF  will  threshold  at  a lower  SNR  than  when  the  noise 
is  uncorrelated.  For  example,  for  y = 0.1,  the  ATF  will  threshold 
at  -32  dB. 


*The  reader  is  cautioned  that  the  assumption  of  noise  coherence,  which 
is  real  and  independent  of  frequency,  is  a highly  artificial  case  and  is 
not  representative  of  real  data.  It  was  chusen  here  primarily  to  indi- 
cate a general  trend  that  could  be  expected  without  becoming  involved 
in  any  particular  noise  model. 
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CASE  III— SINGLE  TARGET  IN  UNCORRELATED  NOISE  PLUS  INTERFERENCE 
The  conditions  for  case  III  are 

6 = 0,  y = 0, 

so  that  equation  (19)-(21)  reduce  to 


<v,(f)  ■ V2(f)  = 1 * gh  + 

(34) 

r lf\  - a-12irftcl  x r - i 2iTf t? 

°w,w2(f)  - e si  ♦ G„  e I . 

Figure  14  shows  the  bias  and  random  errors  as  a function  of  t.  for 
the  following  parameter  values: 

t<--|  = 0 (target  at  broadside) 

tj  varies  form  0 to  1 .6  msec 

Gjj  = 2.5  (interference  is  2.5  times  as  intense  as  the  signal) 

G = 6.32  (SNR  = -8  dB) 
nn 

Also,  for  each  value  of  t, , the  automatic  search  routine  was  initiated 
at  the  true  target's  tinredelay,  that  is,  at  b = 0 

Before  discussing  figure  14  further,  it  should  be  noted  that  the 
resultant  error  function  for  this  case  is  the  sum  of  the  error  function 
due  to  the  true  target  (as  shown  in  figure  7)  and  to  the  interfering 
target  (2.5  times  as  large  as  that  shown  in  figure  7 and  shifted  to 
the  right  by  t.).  In  light  of  the  discussion  concerning  case  I,  the 
ATF,  when  operated  in  its  lowest  mode,  will  seek  the  zero-crossing 
of  the  resultant  error  function  having  the  smallest  bias  error  and  a 
positive  slope. 

The  ATF  in  figure  14  will  settle  on  a time  delay  situated  between 
the  true  target  and  the  interfering  target.  Indeed,  it  may  be  shown 
that  for  small  tj  the  resultant  bias  error  is  approximately  given  by 
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As  t.  increases  further,  this  linear  relationship  will  no  longer 
hold  and^he  bias  error  will  almost  approach  tj.  Then,  in  the  vicin- 
ity of  tj  » 0.8  msec,  the  bias  error  abruptly  changes  sign  due  to 
the  effect  of  the  first  "side  lobe"  of  the  interference  error  function. 
Eventually,  the  bias  error  smoothly  decreases  to  zero  as  tj  approaches 
1.58  msec,  at  which  point  a null  occurs  in  the  error  function  because 
of  the  interfering  target.  It  is  important  to  note  that  the  discon- 
tinuity in  the  bias  error  depends  on  the  search  routine  and  on  the 
initial  bias  error.  If,  for  example,  the  automatic  search  routine 
were  initiated  at  a delay  time  in  the  vicinity  of  tj,  the  ATF  would 
have  continued  to  track  the  interference  and  no  discontinuity  would 
have  occurred.  As  tj  increases  from  1.58  to  2.72  msec,  the  bias 
error  also  exhibits  a cyclic  behavior.  The  maximum  excursions  are 
considerably  smaller,  however,  since  the  2nd  side  lobe  of  the  inter- 
ference error  function  is  smaller.  There  is  no  discontinuity  in  this 
region  because  the  slope  of  the  resultant  error  function  remains 
positive  in  the  vicinity  of  the  origin.  Also  notice  that,  as  for 
case  II,  the  random  error  is  only  a small  fraction  of  the  bias  error. 

CASE  IV--SINGLE  TARGET  WITH  MULTIPATH  IN  UNCORRELATED  NOISE 

The  simple,  but  practical,  multipath  example  considered  for  case  IV 
is  shown  in  figure  15.  In  that  figure,  the  two-element  array  is  situ- 
ated at  a fixed  depth  (D2),  whereas  the  target,  which  is  located  at  a 
fixed  range  (R  and  angle,  G with  respect  to  the  array  normal),  varies 
in  depth  from  200  to  1000  ft  (60.96  to  304.80  m).  The  acoustic  signal 
impinging  upon  the  array  arrives  via  a direct  path  and  a surface- 
reflected  path.  The  problem  is  to  determine  the  bias  and  random  errors 
as  a function  of  the  target  depth  for  various  angles  of  0. 

From  the  geometry,  one  may  readily  compute  the  required  parameters 
listed  on  pages  8 and  9: 

t - t t - L0S1  ' L0S2 

z0  0S1  Z0S2  - C 
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t - ! 

R sine 

1 

si  c 

L0S1 

\ for 

R >>  d, 

t . d 

R sine 

lS2  " C 

L0S2  1 

where 

" o 

. 2 ~ 

1/2 

L0S1 

| + (D1 
L 

~ °2> 

length  of  the  direct  path 

~ ? 

A 9 ' 

1/2 

L0S2  = 

R + (D1 

+ d 2r 

length  of  the  reflected  path 

1 

D-j  = depth  of  the  array 

= depth  of  the  target 

d = sensor  spacing 

Further  particulars  assumed  for  this  case  are 

1 _ * 

6 = - j (negative  due  to  the  180°  phase  shift  at  the  surface 

l0S2  reflection) 

C = 4850  ft/sec  (1478.28  msec) 

d = 3.692  ft  (1 .1253  m) 

D1  = 600  ft  (182.88  m) 

R = 3000  ft  (914.40  m) 

I = 0 (no  interference) 

y = 0 (noise  uncorrelated) 

Gpn  = 6.32  (SNR  = -8  dB) 

The  result  of  simultaneously  solving  equations  (8)  and  (17)  for 
the  bias  and  random  errors  is  shown  in  figure  16.  It  will  be  noted 
that  the  random  error  is  essentially  the  same  for  the  three  angles 
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illustrated.  This  is  merely  a further  demonstration  of  the  fact  that 
equations  (8)  and  (17)  are  only  loosely  coupled  at  small  bias  and  small 
random  errors.  The  bias  error  is,  of  course,  zero  only  when  the  tar- 
get is  at  broadside  (0  = 0)  and  increases  smoothly  as  0 increases. 

The  magnitude  remains  relatively  small,  however,  because  the  difference 
between  and  is  small. 


CONCLUSIONS 

It  has  been  shown  that  the  straightforward  application  of  simple 
control  theory  to  a linearized  model  of  an  analog  BDI  results  in  two 
coupled  transcendental  equations  whose  unknowns  are  the  mean  and  the 
variance  of  the  controlling  time  delay.  The  unique  feature  of  this 
analysis  is  that,  from  the  outset,  this  time  delay  is  treated  as  a 
random  variable  (Gaussian)  and  that  the  BDI  is  linearized  about  an 
arbitrary  operating  point.  This  ensures  that  the  resulting  expres- 
sions remain  valid  at  large  rms  errors  and  allows  investigation  of 
the  ATF  performance  in  the  vicinity  of  the  threshold. 

It  was  also  shown  that  a more  complete  expression  for  the  BDI 
gain  is  one  that  involves  the  total  variation  of  the  average  BDI  out- 
put, that  is,  the  variation  wi  tFi  respect  to  the  mean  and  the  standard 
deviation  of  the  controlling  time  delay. 

The  utility  of  this  analytical  approach  has  been  demonstrated  by 
solving  the  coupled  mean  and  variance  expressions  for  split-beam  ATF 
operating  in  four  adverse  acoustic  environments.  Particular  emphasis 
was  place  on  case  I (Single  Target  in  Uncorrelated  Noise)  since  this 
afforded  a simple  example  of  the  calibration  procedure  of  the  effec- 
tive time  constant  of  the  closed-loop  filter,  as  well  as  of  the  iden- 
tification of  stable  operating  modes. 

This  analysis  has  shown  that  the  estimator  is  unbiased  only  when 
the  ATF  is  operating  in  its  lowest  order  mode,  and  that  the  variation 
of  the  random  error  with  the  SNR  is  identical  to  that  predicted  by 
the  Cramer-Rao  lower  bound  over  most  of  its  operating  range.  At 
-18.6  dB,  the  ATF  exhibits  the  expected  threshold,  whereas  at  an  SNR 
in  excess  of  approximately  +4  dB  the  ATF  exhibits  slightly  higher 
random  errors.  As  the  noise  becomes  partially  correlated,  the  bias 
error  increases  sharply  and  is  considerably  larger  than  the  random 
error.  As  expected,  at  low  SNR,  the  ATF  tracks  the  noise,  which  is 
effectively  located  at  broadside. 
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Cases  III  and  IV  (Single  Target  in  Uncorrelated  Noise  Plus  Inter- 
ference and  Single  Target  with  Multipath  in  Uncorrelated  Noise, 
respectively)  show  the  effect  of  an  interfering  target  and  multipath. 
These  cases  were  included  primarily  to  demonstrate  the  general  use- 
fulness of  the  derived  expressions. 

It  should  be  mentioned  again  that  the  work  in  this  report  was 
confined  to  a two-element  array  in  order  to  focus  on  the  ATF  analysis, 
and  may  readily  be  extended  to  include  any  array  geometry  by  proper 
modification  of  the  equations  in  appendix  A. 

Indeed,  the  computer  program  mentioned  in  the  text  is  actually 
written  for  an  unshaded  50  element  equally  spaced  array.  With  proper 
scaling  to  the  effective  half-beam  centers,  the  primary  effect  of  this 
array  is  to  decrease  the  height  of  the  "side  lobes"  of  the  error  func- 
tion. Consequently,  the  results  for  case  I remain  essentially  un- 
changed (except  for  the  increased  SNR),  whereas  the  effect  of  an 
interfering  target  is  decreased  considerably  from  that  indicated  in 
case  III  even  after  proper  scaling. 
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90'  PHASE  SHIFTER  OR  DIFFERENTIATOR 


Figure  1 . A BDI 
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Figure  3.  Linearized  Form  of  the  BDI  That,  in  the  Steady-State, 
Matches  the  Mean  and  Variance  of  the  True  BDI  Output 


Figure  4(a).  Linear  Model  for  the  ATF  Loop 


Figure  4(b).  Equivalent  Filter  Relating  N'  and  e 
Figure  4.  Model  for  Determining  the  Closed-Loop  Filter 


Diagram  of  the  ATF  Under  Uncorreiated 
Noise  Conditions 


ERRORS  </Juc) 
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Appendix  A 

THE  MEAN  AND  VARIANCE  OF  THE  BDI 

The  mean  and  the  variance  of  the  output  of  the  processor  shown  in 
figure  1 are  derived  here.  This  is  accomplished  by  first  computing 
the  relevant  conditional  expectations  and  then  averaging  them  over  the 
random  time  delay  (e). 

In  figure  1,  the  Fourier  transform  of  the  composite  filter  in 
channel  1 is 


H-j  ( f ) = H(f ) 


(A- 1 ) 


and  in  channel  2 is  (at  a given  e) 


H2(f)  = H(f ) Y(f )e' 


•i  2itfe 


(A-2) 


The  auto-  and  cross-spectral  densities  of  x-j(t)  and  x^t)  may  then  be 
expressed  as 


V/f)  ' |H(f)|2Sx,(f)  <A'31 

Wf)  ■ l¥f>l2Wf)  = |H(f>|2  IVtfli2^)  (A-4) 


2 2 


2 2 


2 *, 


V2(,-e)  ■ Hi(f>Hz(f>  v2(f)  = |H(f)|  v (f)  v2(f) 


i 2-nfe 


(A-5) 


where  G (f)  and  G ,,  (f)  are  the  auto-spectral  densities  and  G (f) 
w1w1  w2w2  w^2 

is  the  cross-spectral  density  of  the  input  sample  functions.  It  should 
be  noted  that  only  the  cross-spectral  density  G (f,e)  depends  on  the 
delay  (e).  X1X2 

For  algebraic  convenience,  the  following  definitions  shall  be  made: 


A-l 
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Vf)  = !H(f)|2  GwlWl(f)  = GXlx/f)  (A-6) 

*2(f)  = |H(f)|2  |Y(f)|2  G (f)  ■ G (f)  (A-7) 

^12(f)  = |H(f)!2  Y*(f)  (f)  = GXlx2(f‘e)e"12lTfe-  (A-8) 

The  conditional  expectation  of  the  output  u(t)  may  now  be  deter- 
mined. At  a given  e. 


U ( t)  = E(u | e)  = X^t)  x 2 ( t ) = R^^CO)  = /GXlx2(f*e)df’  (A_9) 

where  the  fact  that  R (t)  and  G (f)  are  a Fourier  transform 

pair  has  been  used.  The  conditional  expectation  averaged  over  the 
probability  density  (p(e))  of  e will  then  yield  the  average  output 


(mu).  Thus, 

mu  = = y^e  P(£)  E(u|£)  * df  p(0  Gx^x^(f,e).  (A-10) 

Substituting  for  G (f,e)  and  integrating  over  e yields 
xlx2 

mu  = f <l>12(f)  c(2ltf)df  (A"11) 

mu  =/  |H(f)|2  Y*(f)  Gw]w2(f)  c(2*f>df’  (A-12) 

where 

C(2irf ) = J ei27rfe  p(e)de  (A- 13) 


is  the  characteristic  function  of  the  random  variable  (e). 
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The  conditional  expectation  of  u2(t)  may  be  computed  in  a similar 
fashion.  For  this  purpose,  it  is  necessary  to  evaluate  u(t)  u(t  - t). 
Using  the  Gaussian  property  of  the  input  sample  functions,  one  may 
write 


u(t)  u(t  - t)  = P.uu(t)  = x^t)  x2(t)  x] (t  - t)  x2(t  - x) 


(A- 14) 


Ruu(t)  ■ * %iX/ORXjX2(t)  ♦ RXiX2<t)RVi(t>.  (A-15) 

Taking  the  Fourier  transform  yields 

Guu(f-e>  ■ Rx,x2(0)  S(f)  * Gx,x/f>®%x2<f)  + Gx,x2<f-£> 

©G*  (f,e),  (A- 16) 

X1  2 


where  © denotes  convolution. 

The  conditional  expectation  of  u2(t)  may  now  be  obtained  as  follows : 


E(u2|e)  = Ruu(0)  =/Guu(f,e)df 

E(u2|e)  = E2(u  | e)  +y'df[GXlx1(f)  ® Gx  (f) 

+ G (f,e)  © G*  x (f »e )1 

A ^ A ^ Ai A2  j 

since,  from  equation  (A-9),  R (0)  = E(u|e). 


(A- 17) 


(A- 18) 


The  variance  of  u may  be  obtained  from  the  following  relation: 


ou2  = var(u)  = E(u2)  - E2(u) 
ou2  = J dz  p(e)  E(u2|e)  - E2(u). 
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Substituting  equation  (A-18)  into  (A-19)  yields 


*u2  = |y"de  p(e)  E2(u  | e)  - E2(t 


/dfde 

v y (f.e)  © Gv  x (f»eM  * 
T2  12  J 


+ G (f,e)  ® Gv  x (f»e)  * 

Xl  2 12  J 

Finally,  substituting  equations  (A-6)-(A-9)  into  the  above  expression 
and  integrating  over  e yields 

°u2  = Jjf ^fldf2  *i2(fl}  ^12(f2}  {C»(fl  + f2)]  " C(2irV  C(2"V} 

+ /*GNN(f}df,  (A- 20) 


where 


GNN(f)  = ^(f)  © *2(f)  + A <t>l2(x)  4>12(x  - f)  C[2tt(2x  - f)]. 

(A-21) 

A moment  of  reflection  will  reveal  that  the  spectral  density  G^uCf) 
is  the  average  (averaged  over  c)  auto-spectral  density  of  a random 
variable  (N),  which,  at  a given  z,  is  equal  to  (u  - u). 

Equations  ( A- 1 1 ) and  (A-20)  represent  the  desired  expressions  for 
a general  random  time  delay  (e).  Any  additional  simplification  of 
these  results  requires  further  assumptions.  Consider,  for  example, 
that  z is  deterministic  and  has  the  value  e . In  this  limiting 
case,  the  probability  density  of  z assumes  The  form  p(e)  = 6(e  - eQ); 
and  the  characteristic  function  of  e becomes 

C ( n ) = /*d£  6(e  - e )e1ne  = e1neo 

•/  0 • 
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It  is  imnediately  apparent  that  the  first  term  in  equation  ( A- 20 ) 
vanishes.  (Hence,  this  term  arises  due  to  the  randomness  of  e.)  The 
complete  expressions  under  this  condition  are 

mu  - fdf  $12(f)  ei2lTfEo 

ou2  = y*df  [^(O®  $2^  + el2ltfe°  ® $12*(f)  e"i2irfeo]- 

(A- 22) 

Another  interesting  simplification  results  if  the  random  time 
delay  (e)  has  a Gaussian  probability  density  with  a mean  (m)  and 
a variance  (o2).  Then,  the  characteristic  function  becomes 


CCn)  - e,nm 


(A- 23) 


The  term  in  the  brackets  in  equation  (A-20)  may  be  expanded  as 
follows: 

jc j2ir(f-j  + f2)]  - C(2irf1)  C(2uf2)J  = [e  ^ f^2  -l] 

.j^  i2irf1m  -l/2(27rf1a)2]  i2irf2m  -l/2(27rf2a)2] 
j j i2irf1m  -l/2(21If1o)2]  £ i2irf2m  - 1 /2( 2irf 2a ) 2] 


£ (12*^0)"  (i2Trf2a)n 


( A- 24 ) 


-(2ttct;  f.fp 

where  the  series  expansion  for  the  exponential  e has 

been  used.  As  a result,  the  first  term  of  the  variance  expression 
(A-20)  becomes 
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[/d f(12*f)n  *12(f)  ei2irfm  ei2,rfa)2]  . 
n=l  ’ L 


However,  the  term  within  the  brackets  is  nothing  more  than  the  n-th 
derivative  with  respect  to  m of  m , as  the  reader  may  readily  verify 
by  substituting  equation  (A-23)  iHto  ( A- 1 1 ) and  differentiating.  Con- 
sequently, the  complete  expressions  for  a Gaussian  random  delay  (e) 
become 


mu  c /df  *12(f)  el2*fm  e“1/2^27rfa^  (A-26) 

“U2-E#  (^)*AN(f)df-  (A-27> 

n=l 

where 

GNN(f)  - ^(f)  © *2(f)  + / dx  <J>1 2(x)  *12(x  - f)  ei2lTm(2x  - f) 
•e-l/2[2iro(2x  - f)]2. 


(A- 28) 


r 
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Appendix  B 

LINEARIZATION  OF  THE  VARIANCE  OF  THE  BDI  OUTPUT 

The  variance  of  the  BDI  output,  as  given  by  equation  (2)  or  (A-27), 
is  expanded  here  about  some  unspecified  operating  point  at  m and  a . 
For  this  purpose,  a number  of  definitions  and  assumptions  ar£  made.0 
Recall  that  the  mean  of  the  BDI  output  (m  ) is  given  by  (see  A-26) 


„ im  f.*  , t _i2Trfm  -l/2(2irfo)2 

nu(m,a  ) = / df  e e 


Now,  if  one  defines  Kn(m,a  ) as 


K (m,a2)  = [m  (m,a2)l  , 

3m  L J 


(B-l) 


(B-2) 


then  the  variance  of  the  BDI  output  may  be  written  as 

2,  2\  2*2  x 1 4 . 1 6„2 

°u  *m,a  / “ ° a ^2  + "?!  a ^3  + ••• 


+• J GNN(f  ,m,a2)df . 


( B- 3) 


By  numerically  computing  the  higher  order  K's,  as  well  as  GNN 
(f,m,a2),  we  can  see  that  within  the  expected  operating  range 
of  the  ATF  the  infinite  series  is  accurately  represented  by  the  first 
two  terms,  and  GNN  is  relatively  independent  of  m and  o.  The  latter 
is  particularly  true  at  low  SNR's.  Hence,  a good  approximation  of 
equation  (B-3)  is 

ou2(m,c2)  » a2K2  + \ a4  + ^G^f  ,m,a2)df . (B-4) 

Considering  the  variables  to  be  m and  a2,  the  above  relation- 
ship may  be  expanded  about  the  operating  point  at  m and  a to 
yield  0 





TR  5307 


Oy(m,o2)  « a2(mo,o2)  + (m  - mj 


o'  \am  /m0,o0 


t / 2 2, 

+ (o  - oQ) 


L3 (a  )J  mo,ao 


(B-5) 


Using  the  notation  that  Kn0  is  the  n-th  derivative  of  m with  respect 
to  m evaluated  at  the  operating  point  at  m0  and  a0,  we  Hote  that  the 
derivatives  in  equation  (B-5)  may  be  expressed  as 


(— ) 


" 2oo  K10  K20  + ao  K20  K30 


(B-6) 


m .a 
o o 


p°u 

= K2  + 2ojKlo(- 

La(»2] 

J 0 

m rr 

o ' 

m, 

O 0 

22  4 /3^2\ 

+ K20  °o  + ao  K20  (“tJ 

v oO  ' 


O’  0 


V°o 


(B-7) 


where  the  fact  that  Gm^  is  constant  with  respect  to  m and  o has  been 
used.  By  using  equations  (B-l)  and  (B-2),  we  see  that  for  all  m and 


9*S  1 1 

_>  = Ik  and  — ^ = Ik.. 

H<n  23  a{7)  24 


Thereby,  equations  (B-6)  and  B-7)  become 


J 

J 


B-2 
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V°o 


= o^K20(2K10  + oQ  K3Q)  w 2a o K1()K20 


(B-8) 


3(o2) 


= ^10^10  + ao*^30^  + °n  + a °n^in) 


o 20v  20  2 o 40; 


Vao 


* K"  + a 


(B-9) 


where  the  approximations  follow  since  Kx  » a K3  and  K2  » o 1C 
within  the  operating  range.  Substituting  equations  (B-8)  and  (t 
and  af(m  .a*)  into  (B-5)  yields 


u'  0’  o' 


72(m,o2}  « a2K20  + 2 aQK20  + >m*a0)df 


(B-10) 


+ (m  - mQ)  ( 2a0Ki qK20^  + “ ao^K10  + °olW’ 


Upon  canceling  terms  and  making  the  reasonable  assumption  that 

df  >>  \ °oK20  ’ 
we  note  that  equation  (B-10)  reduces  to 

02(m,o  ) ~ a (K-jq  + oqK2q)  + " mo^  2aoK10K20 

+ y*GNN(f,mo‘  ao)df* 


(B-l 1 ) 


B-3 


TR  5307 


TR  5307 


INITIAL  DISTRIBUTION  LIST 


Addressee  No.  of  Copies 

COMSUBDEVGRUTWO  1 
CNO,  0P-02  (2),  -224,  -22,  -942U,  -953,  -981G,  -009T  (2)  9 
CNM,  ASW-14  (2),  -23  3 
DDR&E  (G.  Cann)  2 
NAVSHIPRANDCEN,  CARD  2 
NAVELECSYSCOMHQ,  PME-107  2 
NAVSEASYSCOMHQ,  SEA-06H1 , -06H2,  -660,  -660D,  -660F  (2), 

-660F1  , -660B1 , PMS-393,  -396  10 
DTNSRDC  2 
NELC  1 
NAVUSEACEN  1 
NISC,  Code  20,  2252  (12)  13 
NAVPGSCOL  2 
DDC,  Alexandria  12 
Naval  Coastal  Lab,  Panama  City,  FL  (T.  Watson,  Code  732)  1 


